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Abstract 

In Quantum Mechanics, the so called VT — symmetric bound states ip(r) can live 
on certain nontrivial contours C of complex coordinates r. We construct an exactly 
solvable example of this type where Sturmian bound states exist in the absence of 
any confining potential. Their origin is purely topological: our "quantum knot" 
model uses certain specific, spiral-shaped contours which iV— times encircle a 
logarithmic branch point of tp(r) at r = 0. 

PACS 03.65. Ge, ll.10.Kk, 11.30.Na, 12.90. +b 



1 e-mail: znojil@ujf.cas.cz 



1 Introduction 



For any spherically symmetric single-particle Hamiltonian H = —A + V(\x\ 2 ) acting 
in L 2 (R D ) let us recollect the related radial Schrodinger equation 

~ 4>(r) + ^ t ^ Mr) + V(r 2 ) Mr) = E ^(r) , r e R + (1) 
dr l r l 

with I = (D — 3)/2 + m in the m— th partial wave. At D = 1 one only needs m = 
(for the even parity states) and m — 1 (for the odd parity states) while at all the 
higher dimensions D > 2, the sequence of the angular-momentum indices m in eq. (CQ) 
is infinite, m = 0, 1, 

At any dimension D > 1 one can shift the coordinate r — ► r e = s — i e using 
a constant e > and a new real variable s G -R in eq. ([T|). As far as we know, the 
resulting change of the Hamiltonian H^ radml ^ — > jj(shtfted) j-^g ]3 een g rs ^ proposed 
by Buslaev and Grecchi [1] who choose their parity-symmetric real potential V(r 2 ) in 
a very specific anharmonic-oscillator form. They proved that their manifestly non- 
Hermitian new Hamiltonian }f( sht f ted ) ^ ^g(shifted)\ can De perceived as physical 
since it is strictly isospectral to a one- dimensional double-well oscillator model with 
manifestly Hermitian and safely confining Hamiltonian Q = Q\ Marginally, they 
also noticed (cf. their Remark 4, loc. cit.) that }{( shl f ted ) [ s PT— symmetric, with 
T denoting the complex conjugation operator and with V representing the standard 
operator of parity. Unfortunately, the latter Remark passed virtually unnoticed [2]. 
Only five years later, Bender with coauthors [3l HJ E] proved much more successful 
in re-attracting attention to the VT— symmetric family of the bound-state models 
where, as they conjectured, the spectrum of the energies has a very good chance of 
remaining real and observable. 

In the latter context let us recollect our recent studies [6l[71[8] and, in their spirit, 
replace the Buslaev's and Grecchi's complex integration contour 

C = C {BG) = {r = s - ie | e > , s G R} 

by an element C^ N ' of a family of its generalizations which will be specified below (cf. 
section [2]). In order to simplify the construction (cf. section [3]) and the discussion 
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(cf. sections H]) , we shall only consider the simplified, dynamically trivial model with 
vanishing V(r 2 ), 

-^^ r ) + ^ -^A^ r )= E ^ r )^ reCW. ( 2 ) 

In order to extend the scope of such a model slightly beyond its purely kinematical 
version (cf. also the summary in section [5]), it is easy to add an interaction term 
V(r) = 7/r 2 . Then, the implicit definition of the effective i in (j2J), 

£(£ + l) =1+ (m + £^^ (m + £-ly m = 0,l,... (3) 

enables us to treat our 7 or I = £(7) as a continuous real parameter. Under this 
assumption we intend to show that in spite of absence of any confining force, our 
Schrodinger eq. ([2]) defined along topologically nontrivial integration paths generates 
Sturmian bound states ip(r) £ L 2 (C^ N n at certain discrete sets of 7 and/or £ = £(7). 
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Figure 1: Sample of the curve with N — 1. 

2 Integration paths C^ N > 

In certain complex domains of r, the general solution of any ordinary linear differen- 
tial equation of the second order can be expressed as a superposition of some of its 
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Figure 2: Sample of the curve with N = 2. 
two linearly independent solutions, 

t/j(r) = di) {1 \r) + c 2 ip {2) (r). (4) 

In our present class of models (j2J) let us distinguish between the subdomains of the 
very small, intermediate and very large \r\. In the former case we can choose 

ip^[r) = r e+1 + corrections , ip^ 2 \r) = r~ l + corrections , \r\ <C 1 , (5) 

while we would prefer another option in the latter extreme, with k = \/E in 

■0^' 2 ^(r) = exp (±i«r) + corrections , \r\ 3> 1 ■ (6) 

In between these two regimes, our differential eq. ([2]) is smooth and analytic so that 
we may expect that all its solutions (HI) are also locally analytic. 

In the vicinity of the origin r = our differential equation has a pole. For the 
time being let us assume that our real parameter I in eq. (T5]) is irrational. In such 
a case, both the components of our wave functions (jl]) (as well as their arbitrary 
superpositions) would behave, globally, as multivalued analytic functions defined on 
a certain multisheeted Riemann surface TZ. In the other words, our wave functions 



possess a logarithmic branch point in the origin, i.e. a branch point with an infinite 
number of Riemann sheets connected at this point [9]. 

Separately, one might also study the simplified models with the rational is which 
correspond to the presence of an algebraic branch point at r = which would connect 
a finite number of sheets [9]. Let us briefly mention here just the simplest possible 
scenario of such a type where £(£ + 1) = to that either £ = —1 or £ = 0. In such a 
setting, eq. (0J just separates ip{r) into its even and odd parts so that the Riemann 
surface itself remains trivial, TZ = <D. 

In the generic case of a multisheeted TZ we intend to show that the asymptotically 
free form of our differential eq. (j2J) with the independent solutions ([6]) can generate 
bound states. One must exclude, of course, the contours running, asymptotically, 
along the real line of r since, in such a case, both our independent solutions ■0( 1,2 )(r) 
remain oscillatory and non-localizable. The same exclusion applies to the parallel, 
horizontal lines C^ BG ^ in the complex plane of r. In the search for bound states, both 
the "initial" and "final" asymptotic branches of our integration paths C^ N ' must have 
the specific straight-line form r = ±|s| e ltp with a non-integer ratio tp/n. Thus, we 
may divide the asymptotic part of the complex Riemann surface of r G TZ into the 
sequence of asymptotic sectors 

S = {r = -ige iv \g^>l, cp G (-tt/2, tt/2)} . (7) 

<S±k = {r = — i e ±iklT g e iv \ g 1 , (p E (— 7r/2, tt/2)}, A; = 1,2,.... (8) 

We are now prepared to define the integration contours C^ N \ For the sake of con- 
venience, we shall set all their "left" asymptotic branches C^ le ^ in the same sector 
So and specify r = (s + s ) (1 + is) where s G (— oo, — s ), e > and s > 0. The 
subsequent middle part of C^ N ' must make iV counterclockwise rotations around the 
origin inside TZ while s G (— sq,sq). Finally, the "outcoming" or "right" asymptotic 
branch of our integration contour with s G (so, oo) must lie in another sector 
S2N of TZ, i.e., in the Riemann sheet where the requirement of VT— symmetry [7] 
forces us to set r = (s — s ) (1 — ie). In this spirit, our illustrative Figures [Tj- [3] 
sample the choice of N — 1, iV = 2 and iV = 3, respectively. 
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Figure 3: Sample of the curve with N = 3. 

3 Bound states along nontrivial paths 

It remains for us to impose the asymptotic boundary conditions requiring that our 
wave functions vanish at s — > ±oo. As long as our integration path performs 
N counterclockwise rotations around the origin, this form of the asymptotic bound- 
ary conditions will already guarantee the normalizability of our bound-state wave 
functions ifj(r) G L 2 (<3( N ^ (cf. the similar situation encountered in the models with 
confining potentials [U El [TO]). 

In the subsequent step we may set z = kt and ip(r) = y/z(p(z) in our bound- 
state problem (j2J) with the (by assumption, real) E = n 2 . This reduces eq. (j2J) to 
the Bessel differential equation with the pair of the two well known independent 
special-function (say, Hankel-function [TTJ) solutions which may be inserted in our 
ansatz 

ilj(r)=c 1 ^Hl 1 \Kr) + c 2 V^H{ ) 2 \Kr), v = £+1/2. (9) 
At |argz| < 7i and Rev > —1/2, the asymptotics of its components are given by the 
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respective formulae 8.451.3 and 8.451.4 of ref. [TT] . 

exp 

^«-p[-,(.-^)] +•••). 

This implies that inside the even-subscripted sectors <S2fc our ansatz (Q combines 
the asymptotically growing (and, hence, unphysical) component Hj) 1 \z) with the 
asymptotically vanishing and normalizable, physical component H^\z). Vice versa, 
in all the odd-subscripted sectors S 2 k+i we would have to eliminate, in principle, the 
asymptotically growing Hj^'(z) and to keep the asymptotically vanishing H^(z). 

We may start our discussion of the existence of the localized bound states from 
the straight-line contour C = C^ BG "> = which is all contained in the zeroth sector 
Sq. This immediately implies that with Imr <C —1, the asymptotically vanishing 
solution 

ip^(r) = y/r H^{nr) = const ■ exp(— inr) + corrections 

remains unconstrained at all the real k. Obviously, the spectrum remains non-empty 
and bounded from below. This means that the low-lying states remain stable with 
respect to a random perturbation. A less usual feature of such a model is that its 
energies densely cover all the real half- line M + . This feature is fairly interesting per 
se, although a more detailed analysis of its possible physical consequences lies already 
beyond the scope of our present brief note. 

Our eigenvalue problem becomes not too much more complicated when we turn 
attention to the spiral- or knot-shaped integration contours with N > 0. In such 
a case, fortunately, the exact solvability of our differential equation enables us to re- 
write ansatz in its fully explicit form which remains analytic on all our Riemann 
surface 1Z. Once we choose our "left" asymptotic sector as S , the "left" physical 
boundary condition fixes and determines the acceptable solution on the initial sheet, 

ip(r) =cVr~Hl 2) (Kr), re S . (10) 

After the N counterclockwise turns of our integration path & N > around the origin 
this solution gets transformed in accordance with formula 8.476.7 of ref. [11] which 
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plays a key role also in some other solvable models 

H m = sin(l + m)^ w + ^ ™ ^ 

V / sin 7T7y sin 7T7/ 



Here we have to set m = 2N. This means that the existence of a bound state will 
be guaranteed whenever we satisfy the "right" physical boundary condition, i.e., 
whenever we satisfy the elementary requirement of the absence of the unphysical 
component Hjp(z) in the right-hand side of eq. (!IT|) . 

The latter requirement is equivalent to the doublet of conditions 

2Nv = integer , v ^ integer . (12) 

This means that at any fixed and positive value of the energy E = k 2 and at any 
fixed winding number N = 1, 2, . . ., our present quantum-knot model generates the 
series of the bound states at certain irregular sequence of angular momenta avoiding 
some "forbidden" values, 
M - N 

M = 1,2,3,..., M ^2N,4N,6N,... . (13) 



2N 1 

These bound states exist and have the analytically continued Hankel-function form 
([TO]) if and only if the kinematical input represented by the angular momenta i is 
restricted to the subset represented by formula f[T3"j) . 

Our construction is completed. Once we restrict our attention to the purely 
kinematic model with 7 = 0, we can summarize that at the odd dimensions D = 2p+l 
giving £ = n+p — 3/2 we may choose any index n and verify that formula (1131) can be 
read as a definition of the integer quantity M = (2n+2p— 1) iV which is not forbidden. 
At the even dimensions D = 2p we equally easily verify that the resulting M is always 
forbidden so that our quantum-knot bound states do not exist at V(r) = at all. 

The latter dichotomy appears reminiscent of its well-known non-quantum real- 
space analogue, but the parallel is misleading because in quantum case the freedom 
of employing an additional coupling constant 7 enables us to circumvent the re- 
strictions. Indeed, once we select any dimension D, angular-momentum index m, 
winding number N and any "allowed" integer M, our spectral recipe (1131) may sim- 
ply be re-read as an explicit definition of the knot-supporting value of the coupling 
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constant 

(M\ 2 ( D-2\ 2 
7= v2ivj - { m + —) ■ 

This implies that at non-vanishing 7s, the quantum knots can and do exist at all 
dimensions. 

4 Discussion 

4.1 Complexified coordinates 

In the language of physics, our present construction and solution of a new and fairly 
unusual exactly solvable quantum model of bound states is based on the freedom 
of choosing the knot-shaped, complex contours of integration C. This trick is not 
new [U H] and may be perceived as just a consequence of the admitted loss of the 
observability of the coordinates in PT-symmetric Quantum Mechanics [13J. 

From an experimentalist's point of view, the omission of the standard assump- 
tion that the coordinate "should be" an observable quantity, x G M, is not entirely 
unacceptable since the current use of the concept of quasi-particles paved the way for 
similar constructions. Related Hamiltonians could be called, in certain sense, mani- 
festly non-Hermitian. Still, they are currently finding applications in nuclear physics 
(where they are called quasi-Hermitian [H]). The loss of the reality of the coordi- 
nates is also quite common in field theory where the similar unusual Hamiltonians 
are being rather called CPT-symmetric [15] or crypto- Hermitian |16j . 

In a pragmatic phenomenological setting, the fairly unusual nature of the new 
structures of spectra seems promising. At the same time, the formalism itself is now 
considered fully consistent with the standard postulates of quantum theory. In the 
language of mathematics, the emergence of its innovative features may be understood 
as related to the non-locality in x, i.e., to the replacement of the standard scalar 
product 

(0 I 4>) = / i/j*(x)(j)(x)dx 
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by its generalized, nonlocal modifications [131 HZ| 

(4>\<t>)=J ^*{ x ) ©(^ V) <P(y)dx dy . 

Although this leaves an overall mathematical consistency and physical theoretical 
framework of Quantum Theory virtually unchanged [13] , a new space is being open, 
inter alii, to the topology-based innovations. In principle, they might inspire new 
developments of some of the older successful applications of the formalism ranging 
from innovative supersymmetric constructions [18] to cosmology [19], occasionally 
even leaving the domain of quantum physics [20J . 

4.2 Towards a simplification of mathematics 

The very existence of the present knot-type model is certainly based on nontrivial, 
non-Dirac scalar-product Hilbert-space kernels 6 ^ / (which is our present symbol 
for the quantity denoted as T in ref. [H] or as rj+ in refs. [21] or as a factorized 
symbol CV in ref. [15] or VQ in [22] etc). This being said, a part of the price for an 
enhancement of the scope of the theory lies in the necessity of the explicit study of 
the equivalence transformations between the simple quantum model (which is non- 
Hermitian and VT— symmetric) and its "physical" , Hermitian alternative represen- 
tation^) }j( Hermit ' ian ) . This transformation is not only non-local and complicated 
but also strongly Hamiltonian-dependent. 

One could try to re-establish some parallels with the exceptional Buslaev's and 
Grecchi's oscillators [T] where ff( Hermittan ) happened to be, incidentally, still of a 
local, differential-equation nature. In the similar fortunate cases (cf. also ref. [23] 
for another illustration), the construction of }j( Hermttian ) remained reducible to the 
construction of a certain effective interaction V^ e ^\r) e JR. 

Partially, a similar technique of simplification can be also applied to our present 
model with V = 0. Indeed, although it seems that all the variations of the spectrum 
must be attributed to the mere kinematical aspects of the model in question, one 
could employ the philosophy of ref. [B] to imagine that the nonequivalent effective 
interactions V^ e ^^(r) could be also attributed to the nonequivalent knot-shaped in- 

10 



tegration contours C^ N \ An implementation of such a perspective would rely upon 
the conformal mapping 

r = -ie ix/2 . (14) 

It simply maps each Riemann sheet in r on a stripe in complex x— plane where 3?e x 
plays the role of the angular coordinate while 3m x simulates the radial coordinate. 
In the other words, the polar-coordinates-resembling mapping (fT4"l) can be used to 
compress several sheets of a multisheeted Riemann surface 1Z onto the single reference 
complex plane of the auxiliary complex variable x. In this manner an effective 
potential would emerge quite naturally (cf. also ref. [12] in this respect). 

5 Outlook 

We can summarize that our present example belongs to such a type of quantum 
theory where 

• "coordinates" r may complexify and lose their immediate observability, 

• the topological structure of paths C can often be simplified via suitable confor- 
mal changes of variables [6j [7J [8] , 

• a "suitable" [TJJ [15j [2TJ [21] redefinition of the inner product in the Hilbert 

space is needed, 

• complex potentials V(r) G(T are allowed. 

Let us add, finally, a few remarks on the latter, nontrivial and appealing point. 
Firstly, one must imagine that our present, exactly solvable N > 1 quantum- knot 
bound-state problem would become purely numerical after its immersion in virtually 
any external confining potential. Secondly, the survival of the reality of the new 
bound-state spectra (i.e., of the measurability of the energies) would have to be 
individually secured or proved as the necessary condition of the physical acceptability 
of the models with V(r) ^ 0. 
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The latter point has already been properly emphasized by Bender and Boettcher 
[I] who choose the one-parametric family of the one-dimensional power-law forces 
y( BB )(x) = x 2 (ix) s and supported the conjecture E E IR by the persuasive numer- 
ical and semiclassical arguments. The more rigorous (and unexpectedly difficult) 
proof only followed several years later [25J. It is worth noticing that the latter proof 
found also its most natural formulation in the general D— dimensional setting with 
D > 1 and/or with the variable real t. 

One of the merits of working with the power-law-dominated confining potentials 
was that the asymptotic behavior of the related wave functions can be analyzed 
easily (cf. also refs. [U, 0] and [10] where such an analysis has been performed and 
discussed in detail). One of the less expected results of these studies (revealed, tested 
and verified also by several other groups of authors [5l[T2j[26] during 1998 - 1999) has 
been the observation of the explicit contour-dependence of the spectra. In fact, this 
discovery, so important in our present context, is even older since its first indications 
can even be traced back to the pioneering 1993 letter by Bender and Turbiner [27] . 
In our papers [61 [7] we emphasized and extended this observation by contemplating 
the "usual" potentials in combination with the "unusual", topologically nontrivial, 
spiral-shaped contours C = C^ N ' of complex coordinates. 

On this background, the exact solvability of the present, force- free example came 
as a surprise. It could be perceived as an important new illustration of the changes 
of the energy spectra mediated by the use of the non-equivalent integration contours. 
Still, in the nearest future, a more numerically-oriented return to the more realistic 
problem of the contour-shape-dependence of the spectra generated by some nontrivial 
VT— symmetric potentials in eq. seems unavoidable. 
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Figure captions 

Figure 1. Sample of the 
Figure 2. Sample of the 
Figure 3. Sample of the 



curve with N = 1. 
curve C (iV) with N = 2. 
curve with JV = 3. 
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